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Abstract. Recent studies of photon-generation mechanisms in the color-superconducting Color-Flavor Locked 
(CFL) phase of dense quark matter have found 7 -ray emissivities in excess of ~ 10®° erg cm“° s“^ for temperatures 
in the 10-30 MeV range. We suggest that this property can trigger 7 -ray bursts (GRBs) and associated fireballs 
at the surface of hypothetical hot (newly born) quark stars with an energy release of up to 10 ^® — 10 ®° erg within a 
fraction of a millisecond. If surrounded by an accretion disk following its formation, the star’s bursting activity can 
last from tens of milliseconds to hundreds of seconds releasing up to 10®^ erg in total energy. We discuss typical 
features of observed GRBs within our model and explain how quark stars in the CFL phase might constitute 
natural candidates for corresponding inner engines. 
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1. Introduction 

Nuclear matter at high density and small tempera¬ 
ture (T) is expected to exhibit color-superconductivity 
(CSC), induced by quark pairing and condensation at 
the Fermi surface, with energy gaps A~100MeV (Rapp, 
Schafer, Shuryak, & Velkovsky 1998; Alford, Rajagopal, 
& Wilczek 1998) and associated critical temperatures 
Tc—0.6A, above which thermal fluctuations destroy the 
condensate (for a review, see, e.g., Rajagopal & Wilczek 
2000). The experimental relevance of such matter mostly 
pertains to astrophysical objects, in particular compact 
stars. Besides its impact on the equation of state (Lugones 
& Horvath 2002, Alford & Reddy 2003), CSC could af¬ 
fect the emission spectrum of compact stars as encoded in 
its electroweak properties (Jaikumar, Prakash, & Schafer 
2002; Reddy, Sadzikowski, & Tachibana 2003; Vogt, Rapp, 
& Ouyed 2004). If CSC extends to the surface of a hy¬ 
pothetical quark star, pulsed photon emission has been 
suggested as a mechanism for a fireball in the 7 -ray burst 
(CRB) context for the 2-flavour superconductor (2SC) in 
(Ouyed & Sannino, 2002)^. Whereas at moderate densities 

Send offprint requests to: ouyed@phas.ucalgary.ca 
^ In the 2SC phase, up and down quarks pair into a color- 
antitriplet, leaving the quarks of the remaining color unpaired, 
and with five of the eight gluons acquiring a mass. The 3 mass¬ 
less gluons possibly bind into light glueballs subject to fast 
decays into photons (Ouyed & Sannino 2001). 


(and T=0) the existence of the 2SC is still under debate 
(Alford & Rajagopal, 2002; Buballa et al. 2004) - it could 
be superseeded by, e.g., crystalline phases (Alford, Bowers, 
& Rajagopal, 2001; Rapp, Shuryak, & Zahed, 2001) -, the 
Color-Flavor Locked (CFL) phase is the favored ground 
state at sufficiently high density (Alford, Rajagopal & 
Wilczek, 1999). In CFL matter all three quark flavours 
(up, down and strange) have a mass which is negligi¬ 
ble compared to the (quark) chemical potential, /r,, so 
that they participate equally in the color condensation, 
breaking the full (local) color-symmetry and one (global) 
SU{A) chiral symmetry. Therefore, the zero-temperature 
CFL phase is electrically neutral without electrons, col¬ 
orless, and its low-lying excitations are characterised by 
8-1-1 Goldstone bosons (due to chiral and baryon-number 
symmetry breaking). 

In previous work (Vogt, Rapp, & Ouyed, 2004), here¬ 
after VRO, we have explored photon emission and ab¬ 
sorption mechanisms in the CFL phase. Based on the 
Goldstone boson excitations (“generalised pions”) of the 
broken chiral symmetry we have employed a Hidden Local 
Symmetry formalism including vector mesons (“gener¬ 
alised p-mesons”) to assess photon production rates and 
mean free paths at temperatures suitable for CFL. We 
have found that above T ~ 5-10 MeV, the emissivities 
from pion annihilation, tt+tt" —> 77 and tt+tt" —> 7 
(rendered possible due to an in-medium pion dispersion 
relation), dominate over those from conventional electro- 




2 


R. Ouyed et al.: Fireballs from Quark Stars in the CFL Phase 


magnetic annihilation, e“'"e“ —> 77 . Given the very small 
photon mean free path at the temperatures of interest, 
r=10-30 MeV (see Fig. 3 in VRO), emission and absorp¬ 
tion are in equilibrium and reflect the surface temperature 
of the star with a flux corresponding to that of a blackbody 
emitter. Therefore the pertinent cooling will largely de¬ 
pend on how efficiently heat conducted by the Goldstone- 
bosons from the interior can be emitted from the surface 
in form of the thermalized photons. The objective of the 
present article is to (i) investigate the astrophysical con¬ 
sequences of such photon production mechanisms^ for the 
early cooling history of hot CFL stars (Sect. |21) and (ii) 
evaluate whether the resulting fireballs can drive/power 
GRBs (Sect. El. Conclusions are given in Sect.0 


2. Early Cooling of CFL stars 

2.1. Plasma photon attenuation 

The surface emissivity of photons with energies below 
hujp ~23 MeV (wp: electromagnetic plasma frequency) is 
strongly suppressed (Alcock, Farhi, & Olinto 1986; Chmaj, 
Haensel, & Slomihsli 1991; Usov 2001). As shown in VRO, 
average photon energies in CFL matter at temperature T 
are ~3T. Therefore, as soon as the surface temperature 
of the star cools below Ta = hujp/3 ~ 7.7 MeV, we con¬ 
sider the photon emissivity to be shut off (“attenuated”), 
i.e., photon emission only lasts as long as the star cools 
from its initial temperature Tq to Ta. With this in mind 
we proceed to study the thermal evolution of the star. 


2.2. Heat transfer and thermal evolution 

We describe the thermal evolution by a diffusion equation, 


dt dr \ dr) ' 


( 1 ) 


where Cv and k are the specific heat and thermal conduc¬ 
tivity of the star matter, respectively. In the CFL phase, 
these are dominated by massless Goldstone bosons, as 
evaluated in Jaikumar et al. (2002), 


= 7,8 X lO^® T^erg cm'^ , (2) 

0 

and in Shovkovy & Ellis (2002), 

K = 1.2 X 10^^ ^Mev ^GB erg cm“^ s“^ K~^ (3) 

^ We note that our study is different from those involving the 
e''"e“ emission above the surface of CFL stars (see, e.g.. Page & 
Usov 2002 and references therein). The latter scenarios reside 
on the fact that the color-supercondnctive matter of the star 
carries a charge so that an e“ abundance can build up thereby 
generating a critical electric field that in turn produces and 
emits e^e“ pairs and photons. This is quite different from our 
model which is based on photon generation inside the star with 
negligible (e^e“) emission. 


with a Goldstone-boson mean free path (in [cm] in Eq. ®) 


^gb{T) = 


4(21 -8 In 2) 
ISv^ TT TmsV 


exp 


\\Z1 MeV J 


cm. 


(4) 


Effects of a crust have been ignored in the present study 
of the star’s cooling history (see discussion in § 1231 ). 

In dimensionless units, the diffusion equation becomes 


df^ a d (~ 

j 


(5) 


with f = r/R,T = T/Tq,t = t/{R/c), Xqb = Xqb/R and 
a numerical coefficient a=0.512 (c is the speed of light). 
As initial condition we assume the temperature of the star 
at T = 0 be uniform, T{t = 0) — Tq- We further have to 
impose boundary conditions for the temperature gradient 
in terms of the heat flux, F{f) = —XoBdT/df, at the 
center and the surface of the star, F{f = 0) = 0 and 
F{f = 1) = T"*. Due to the large opacity we expect that 
the photons to be thermalised upon leaving the star. Thus, 
we assume black-body radiation to be a good approxima¬ 
tion (supported by VRO). The energy flux per unit time 
through a spherical shell at a radius f is evolved by taking 
into account the gradient of temperature on both sides of 
the shell. For the outermost shell, the flux exits the star 
via thermal photons, implying that the cooling curves re¬ 
sult from photon emission only. 

Numerical solutions to Eq. ®, using a finite-difference 
method, are summarized in Eigs.^andlSlfor a quark star 
of radius i?=10 km. Fig.^shows the temperature vs. dis¬ 
tance from the star’s surface, R—r, indicating that T drops 
to Ta in less than 0.1 ms over a shell of thickness less than 
~1 km. Neglecting heat transport, the cooling timescale 
may be roughly estimated from Cy{dTjdt)=—r=—L j/XV 
{/XV=^ttR^/XR-. volume of the cooling shell of thickness 
/XR). Assuming a blackbody luminosity, L = 
one finds tcooi — 0.1ms x Ai?km x (AT/T), (AT = 
To — Ta), approximately reflecting the numerical results. 
Note that our estimate improves toward small T and large 
A for which Xqb is larger, so that neglecting the delay 
time due to heat transport (Eq. Q) is better justified; 
e.g., for To=20 MeV and A=150 MeV (in which case 
AGB = 10~^kin, reaching 1 km at T~15 MeV), Ai?km—1 
(lower panel in Fig. with (AT/T)~1 one recovers 
Cooi— 0.1 ms, consistent with the lower panel of Fig. 13 

For identical Tq, a larger gap implies cooling deeper 
into the star, but slower in terms of the reduction in sur¬ 
face temperature, cf. Fig. |3 Again, this follows from the 
larger Agb, for which heat emerges from deeper in the star 
and thus provides a larger energy reservoir. 

When increasing the initial temperature, Tq, from 
20 MeV to 30 MeV, the surface cooling does not change 
radically. The main difference is that the cooling of the 
bulk sets in earlier so that the temperature gradient be¬ 
tween surface and bulk is washed out in shorter time. 

We recall that contributions from neutrino {v) emis¬ 
sion are not included in our analysis. By assuming uniform 
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Fig. 1. Radial temperature profiles of CFL stars of ra¬ 
dius R according to solutions of the diffusion Eq. o 
(r: distance to the star’s center) for two different gaps, 
A=50 MeV (solid lines) 150 MeV (dotted lines), and ini¬ 
tial temperatures To=30 MeV (upper panel) and 20 MeV 
(lower panel). In the lower panel, the A=50 MeV case 
cools to Ta=7.7 MeV (lower limit of the plots) in less than 
40 /rs. 


star density and temperature as our initial state the con¬ 
sequences as shown in VRO are (i) v emission from the 
surface is usually negligible in comparison to photon emis¬ 
sion in the temperature region we are considering (see §3.4 
in VRO for more details), and (ii) v emission from the bulk 
will set in once the neutrino mean-free-path becomes com¬ 
parable to the star radius which is expected to occur for 
temperatures below ^5 MeV (Reddy et al. 2002); above, 
the neutrinos are essentially trapped inside the star. This 
point is further discussed in the conclusion. 

We also note that photons will be redshifted as they 
stream outwards due to the star’s gravitational potential. 
First order General Relativity effects on the emissivities 
can be introduced through redshift factors expressed in 
terms of the star mass and radius. We estimate the perti¬ 
nent redshift reduction to be of the order of 10% (cf. Frolov 
& Lee 2005 for a recent analysis). While the gravitational 
redshift will degrade the total energy emitted it should 
not affect the qualitative features of our results. 

3. Application to Gamma Ray Bursts 

The total energy released during cooling is computed from 
Fig. 121 as Etat = / L~fdt. Thus, a CFL star with initial 
temperature of 10-30 MeV can release on average 10^®- 




Time (ms) 


Fig. 2. Time dependence of surface temperature (lower 
panel) and photon luminosity (upper panel) of a CFL star 
for two initial temperatures, ro=20 and 30 MeV (lower 
and upper lines, respectively). Solid (dotted) lines are for 
A=50(150) MeV. 

1050 g]-g -vvithin ~0.1 ms as it cools to Ta. We have at 
hand an engine that could be driving observed GRBs. 

3.1. Accretion disk and temporal variability 

Our main idea for the following is that the input energy to 
the engine is provided by infalling matter from an accre¬ 
tion disk. The energy released by the hot CFL star trans¬ 
lates into (0.01 — \)Mq(? s“^ in accretion energy. These 
values are reminiscent of hyperaccreting disks (Popham, 
Woosley, & Fryer 1999) suggesting that if the latter is 
indeed formed around a CFL star it could reheat the sur¬ 
face of the star via accretion keeping the engine active 
for a much longer time and providing higher energies as 
compared to collapse events only. More precisely, the total 
available energy is related to the disk mass, Mdisk, by 

Et — jyMdiskC^ H- ^^AEb , (6 ) 

Tor 

where the first term accounts for the release in gravita¬ 
tional binding energy of the accreted matter, with ry^O.l 
(e.g., Frank, King, & Raine 1992). The second term rep¬ 
resents the binding energy released in the conversion from 
nuclear to CFL matter (AEb per nucleon). It is presum¬ 
ably emitted via neutrinos (Glendenning 1997) so that 
Et— ryMdiskC^ remains available to photon production. 

We are modelling the effects of the accreted material 
on the CFL star in a very simplistic approach, as a random 
increase of the CFL surface temperature Tg (heating) to a 
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Fig. 3. Model results for the effects of a hyperaccreting disk on the star luminosity; shown is the 7 -ray burst emission 
versus time for three different time resolutions (64 ms, 256 ms, and 512 ms from left to right) representative of GRBs 
detectors (www.batse.msfc.nasa.gov/batse/grb/lightcurve/). The upper (lower) panels are for a total accreted energy 
Et=Q.001Mqc^ (i?T=0.01M©c^) with T„iax=30 MeV. 
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Fig. 4. Same as in Fig. O but for Ti„ax=20 MeV. More sample of lightcurves can be found at 
http: // WWW. capca. ucalgary. ca/ 


value Tpeak with a uniform distribution between the lower 
limit, Ta—7.7 MeV, and an upper limit, Tmaa:=15-30MeV. 
This roughly corresponds to expected average disk tem¬ 
peratures up to Tiiisk—75 MeV or so (Popham et al. 1999); 
we neglect cooling processes during accretion. Note that 
for accreted material with temperature below Ta there is 
no bursting and the engine remains shut off until further 
accretion drives Tg above Pa- In these cases the accretion 
proceeds for longer than the average accretion timescale. 


For the time intervals of accretion, Ataccr, we assume 
the free-fall time for a hyperaccreting disk as the relevant 
scale, iff= 1 / V47rGpdisk —1 ms, with pdisk— 10 ^^ g cm“^ the 
accretion disk density and G the gravitational constant. 
The accretion time intervals are sampled stochastically 
between 1 and 10 ms. The time scale, At.y, of the sub¬ 
sequent 7 -ray burst (cooling) follows from the diffusion 
Eq. @ as the time it takes to cool the surface to Ta- 
During this time the accretion energy is transformed to 
photon energy which in a simple form can be written as 
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Fig. 5. 7 -ray burst emission spectra versus time for a CFL star with accretion disk of mass O.OIMq (or Et = 
O.OOIMqc^; upper panel) and O.IMq (or Et = O.OIMqc^; lower panel), with Tniax=15 MeV and time resolution 
1024 ms. More examples can be found at http://www.capca.ucalgary.ca/ 


L..yAt.y=? 7 M(jiskC^Ataccr- Thus, each event (reheating and 
cooling) lasts for Atevent=Ataccr + for given A and 
Tpeak (we set A=50 MeV). Consequently, each episode 
lasts for a few milliseconds consisting of a linear increase 
of Ts to Tpeak followed by a rather sharp decrease to Ta. 

For a given disk mass (see S id.511 . the simulation is car¬ 
ried out until the total available energy, i?T—O.lMdiskC^ 
(Eq. ®), is consumed. The resulting variability (dis¬ 
played in Figs.|2|and0]for i?T=O.OOlM 0 C^, O.OIMqc^ and 
typical time resolutions of CRB detectors) share interest¬ 
ing similarities to those observed in CRB data, see, e.g., 
http://www.batse.msfc.nasa.gov/batse/grb/lightcurve 
Note that our spectra correspond to the activity directly 
at the engine while the CRB curves represent emission 
beyond the transparency radius (see, e.g., Piran 2000). 
Nevertheless, the engine activity should reflect the 
variability at large distances. 

3.2. Total duration 

Our simulations of the accretion suggest engine activities 
varying from milliseconds to hundreds of seconds depend¬ 
ing on Mdisk and Tmax (Figs. 0 and 0. The average du¬ 
ration can be estimated as tdur—Ataccr x {Et/E.^) with 
average burst energy and accretion timescale erg 

and Ataccr —5 ms, respectively. One obtains tdur—10 s for 
a CFL star surrounded by a O.IMq disk. For this disk mass 
(corresponding to Ft^O.OIMqc^) the engine activity can 
last for hundreds of seconds if Tnax is sufficiently low (e.g., 
15 MeV as in Fig.^J. In general, for a given Mdisk, smaller 


Tmax imply smaller average energies released per episode 
(i.e. smaller accretion rates), and thus longer lifetimes. 

We also note that the variability observed in the first 
few seconds of the Ft= 0 . 01 MqC^ simulations is very sim¬ 
ilar to the Ft^O.OOIMqc^ case which lasts for ~2 s, recall 
Figs.^landQI Similarly, the variability seen in the first few 
seconds of a i?T=O.lM 0 C^ simulation (not shown here) 
bears many similarities to the Et=0.Q1Mq(/ case. This 
“self-similar” behavior is due to the duration being con¬ 
trolled by accretion timescales Ataccr>At.y, in connection 
with the stochastic accretion process whose main effect is 
to spread the events (peaks) randomly in time. 

3.3. Total energy 

In our model up to 10% of MdiskC^ can be transformed 
into fireball energy. E.g., if an accretion disk with mass 
O.IM 0 surrounds a CFL star following its formation (see 
S 13.51) . up to 10 ^^ ergs of energy can be released in 7 -rays. 
Larger disks lead to higher energy release but in this case 
the star is likely to become a black hole in the process. 

3.4. Baryon loading and beaming 

The intense, localised (source size <100 km) and brief ex¬ 
plosion implied by the observed GRB fluxes imply the 
formation of a e“'"e“-photon fireball. Furthermore, most 
of the spectral energy in GRBs is observed at >0.5 MeV, 
so the optical depth for 77 —> e+e" processes is very large. 
Any photon generated above 0.511 MeV will be degraded 
to below 0.511 MeV via the 77 ^ e+e” process. This 
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is the so-called compactness problem leading to a ther- 
malised fireball (Ruderman 1975) unlike the optically thin 
spectra observed in GRBs. 

The compactness problem can be resolved if the emit¬ 
ting matter is moving relativistically toward the observer. 
In this case the relative angle at which the photons col¬ 
lide must be less than the inverse of the bulk Lorentz 
factor, 1/r, to effectively reduce the threshold energy for 
e+e” pair production (Goodman 1986). One can show 
that r > 100 is required to overcome the compactness 
problem (Shemi & Piran 1990; Paczyhski 1990). This rel¬ 
ativistic outflow presumably arises from an initial energy 
Eq imparted to a mass Mq E^jt? close to the central 
engine (cf. Meszaros 2002 for more details). 

Assuming the star is bare^ in our model the 
emitted photons are likely to interact with particles 
from the accreting material. To enable acceleration to 
Lorentz factors above 100 the photon energy must 
be imparted on particles representing on average less 
than 0.1% of the total matter accreted per episode, 
since r=7?macc.C^Ataccr/WejecC^=?7(Waccr/TOejec)- In Other 
words, we require most of the infall material to convert to 
GFL matter with only a small portion of it to be ejected. 
This seems not unreasonable as most of the accreted parti¬ 
cles will instantly deconfine to quark matter upon contact 
with the star’s surface (see, e.g., Weber 2005). 

The amount of mass ejected will vary from one episode 
to another leading to a spread in the Lorentz-factor dis¬ 
tribution as already indicated in Figs. and A fast 
loaded fireball that was injected after a slower one will 
eventually catch up and collide, converting some of the 
kinetic energy of the shells to thermal energy. This is remi¬ 
niscent of the popular internal-shock model for GRBs (see, 
e.g., Kobayashi et al. 1997; Piran 2005) where a succession 
of relativistic shells collide and release the shock energy via 
synchrotron emission and inverse Compton scattering. 

The intricate details of the accretion-ejection cycles are 
beyond the scope of this paper. Qualitatively, upon tak¬ 
ing into account the star-disk magnetic field, we expect 
accretion to be channelled toward the polar regions as il¬ 
lustrated in Fig. El A refined scenario could further include 
hot bursting spots on the surface of the CFL star emit¬ 
ting 7 -rays; the outgoing photons interact with the parti¬ 
cles from the accreting matter ejecting part of it. Thus, in 
our model, the loaded fireballs emanate mostly from the 
polar regions with collimation reinforced by the magnetic 
field (e.g., Fendt&Ouyed 2004). When a hot spot cools 
below Ta, further accretion is triggered leading to another 
photon-bursting spot on the star as pictured in Fig. El The 

® It has been suggested that the plausible surface depletion 
of s quarks induces a Coulomb barrier allowing the presence of 
a crust (Usov 2004). Nevertheless a crust should be blown away 
by the extreme radiation pressure induced by the extreme tem¬ 
perature (tens of MeV) of the star immediately following its 
birth and during its early history (e.g. Woosley&Baron 1992). 
We thus assume that the star is bare for most of the burst¬ 
ing era where the high temperatures are maintained by the 
infalling disk material. 



Fig. 6. Illustration of the intermittent accretion-ejection 
mechanism in our model. The upper panel shows the disk 
material channelled by the magnetic field toward the po¬ 
lar regions reheating the surface temperature to Tpeak (the 
flashy spots). When Tpeak > Ta, a photon burst is trig¬ 
gered ejecting particles from the accreting material. From 
a viewpoint along the polar axis (depicted by a cross in the 
lower panel), the multiple flashes on the surface of the star 
illustrate the locations of previous bursts (or ejections) in¬ 
duced by previous accretion episodes. The spatial spread 
of these flashes portrays the randomness in the accretion 
and cooling timescales. 

time delay between two successive spots is related to the 
randomness in the accretion (reheating) and the cooling 
timescales as described in § EH This process continues 
until Tt is consumed, which in some cases can support 
up to hundreds of episodes (subjets), as in Fig. El 

3.5. Formation scenarios/sites 

In scenarios for the formation of CFL stars massive pro¬ 
genitors are naturally favoured since they are more likely 
to lead to compact remnants with high enough densities in 
the core for a (phase) transition to quark matter to occur. 
This is in line with afterglow observations which have pro¬ 
vided several hints that some GRBs are associated with 
massive star progenitors (Galama et al. 1998, van Paradijs 
et al. 1999). If the quark star forms immediately follow¬ 
ing a supernova (SN) explosion of a massive progenitor, 
a thick disk (Mdisk > O.IMq) is expected to form from 
fail-back material. The faster fireballs from the CFL star 
will catch up with the preceding supernova shell as to en- 
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ergise it. In fact, the energetic subjets from the CFL star 
are capable of destroying the original symmetry of the 
expanding SN shell. Another possibility involving massive 
progenitors is that of a collapsar-like event (Woosley 1993) 
where a quark star is formed instead of a black hole. In 
this scenario as well one would favour thick and massive 
disks to form around the star. 

The hadron-quark transition may also happen long af¬ 
ter the supernova explosion as in the Quark-Nova model 
where a neutron star (originally formed from a massive 
progenitor) can reach quark-matter densities through ac¬ 
cretion from a companion, or due to spin-down (Ouyed, 
Dey, & Dey 2002; Ouyed et al. 2004). In some cases, the 
phase transition may take tens of millions of year to oc¬ 
cur. During the explosion up to O.OIMq of material can be 
ejected (Keranen, Ouyed, & Jaikumar 2005) which later 
turns into a hyperaccreting disk (Keranen&Ouyed 2003). 

3.6. Bimodality: Short and long duration GRBs 

The durations of GRBs observed by BATSE show a bi- 
modal distribution, which has led to a classihcation of 
GRBs into short (with tgo < 2 s, where tgo is the decay 
time to 10% intensity) and long (tgo > 2 s) (Kouveliotou et 
al. 1993; McBreen et al. 1994). It has been suggested that 
different underlying engines are operative, the short ones 
being related to binary neutron star mergers, and the long 
ones to the collapse of massive stars (see Meszaros 2002, 
and references therein). However, recent comparisons of 
short GRBs lightcurves to the first few seconds in long 
GRBs indicate that the two classes could be similar (e.g., 
Nakar&Piran 2002; Ghirlanda et al. 2004). 

As noted in S lim our model provides a kind of “self¬ 
similar” behaviour where the variability in the first few 
seconds of a long duration GRBs (thick disk) is remi¬ 
niscent of that of a short GRB (thin disk), both driven 
by the same engine. As for the observed bimodality, it is 
possible that the formation scenarios discussed in § im 
(namely, collapsar-like events and Quark-Nova explosion) 
lead to 2 distinct classes of accretion disk (thick with 
Afdisk > O.OIMq and thin with Mdisk < O.OIMq) sur¬ 
rounding the GFL star. 

Interestingly, it has been argued recently in the litera¬ 
ture that the bimodality in GRB duration originates from 
discrete emission regions (subjets) in the GRB jet (Toma, 
Yamazaki, & Nakamura 2005 and references therein). In 
this model the multiplicity of the subjets (rig) along a line- 
of-sight differentiates between short GRBs {ug 1) and 
long GRBs {rig ^ 1). Many subjets have to be randomly 
launched by the central engine for the model to work. 
This bears a close resemblance to the picture we devel¬ 
oped here (cf. Fig.j^l and warrants more detailed studies 
of accretion-ejection within our model. 

4. Conclusion 

We have studied consequences of previously calculated 
photon ( 7 -ray) emission mechanisms in GFL matter for 


the early cooling of GFL stars. Based on the notion that 
pertinent emissivities essentially saturate the blackbody 
limit for temperatures T ~ 10-30 MeV (due to underlying 
processes involving the Goldstone bosons of GFL matter), 
we have solved a diffusion equation and found that each 
photon burst can release an average energy of ~ 10 ^® erg 
during a fraction of a millisecond. We have suggested a 
schematic model within which the time variability and the 
long activity (up to hundreds of seconds) of a GRB engine 
is driven by a surrounding hyperaccreting disk resulting 
from the formation process of the GFL star. This model 
reproduces several features of observed GRB spectra. 

In our simplified picture we assume that photons and 
neutrinos are emitted at the same temperature (see sec¬ 
tion 3.4 in Vogt et al. for a discussion). This certainly 
is far from a realistic situation where temperature gradi¬ 
ents will be present and where the neutrinosphere will be 
buried deeper in the star at hotter temperatures than the 
photosphere. In this case neutrinos will dominate the cool¬ 
ing leading to timescale shorter than the 0.1 ms range we 
found when only photons are at play. This is of particular 
importance in the early stages immediately following the 
formation of the star. On the longer timescales, and within 
the GRB context, the reheating mechanism (induced by 
the infalling accretion disk material) affects mostly the 
surface layers and thus specifically the photosphere re¬ 
gion. A decoupling between the photosphere and the neu¬ 
trinosphere is to be expected justifying the neglect of neu¬ 
trinos cooling in the subsequent photon bursts/episodes. 

For a more quantitative description of the complex en¬ 
ergetics and dynamics involved in our model, advanced 
numerical simulations should be performed accounting for 
general relativistic effects, as well as the star-disk mag¬ 
netic field. Such studies will be similar to what has been 
done in the case of black-hole-accretion-disk systems (De 
Villiers, Staff, & Ouyed 2005), with the black hole being 
replaced by a GFL star. 
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